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We have generalized the BOUND and MOLSCAT packages to allow calculations in basis sets 
where the monomer Hamiltonians are off-diagonal and used the new capability to carry out bound- 
state and scattering calculations on 3 He-NH and 4 He-NH as a function of magnetic field. Following 
the bound-state energies to the point where they cross thresholds gives very precise predictions of 
the magnetic fields at which zero-energy Feshbach resonances occur. We have used this to locate 
and characterize two very narrow Feshbach resonances in 3 He-NH. Such resonances can be used to 
tune elastic and inelastic collision cross sections, and sweeping the magnetic field across them will 
allow a form of quantum control in which separated atoms and molecules are associated to form 
complexes. For the first resonance, where only elastic scattering is possible, the scattering length 
shows a pole as a function of magnetic field and there is a very large peak in the elastic cross section. 
For the second resonance, however, inelastic scattering is also possible. In this case the pole in the 
scattering length is dramatically suppressed and the cross sections show relatively small peaks. The 
peak suppression is expected to be even larger in systems with stronger inelasticity. The results 
suggest that calculations on ultracold molecular inelastic collisions may be much less sensitive to 
details of the potential energy surface than has been believed. 
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I. INTRODUCTION 

Over the last 5 years, it has become possible to con- 
trol the behavior of ultracold atomic gases by tuning the 
interactions between atoms using applied magnetic fields 
Q, Q . Notable successes have included the controlled im- 
plosion of Bose-Einstein condensates (jl and the produc- 
tion of molecules in both bosonic 0,|Elail3 an d fermionic 
|E El [HI El quantum gases. Long-lived molecular Bose- 
Einstein condensates of fermion dimers have been pro- 
duced IT^. Il3l Il4| , and the first signatures of ultracold 
triatomic |15| and tetraatomic [HI molecules have been 
observed. The new capabilities in atomic physics have 
had important applications in other areas: for example, 
the tunability of atomic interactions has allowed explo- 
ration of the crossover between Bose-Einstein condensa- 
tion (BEC) and Bardeen- Coop er-Schrieffer (BCS) behav- 
ior in dilute gases (II 111 ITif . 

In parallel with the work on atomic gases, there have 
been intense efforts to cool molecules directly from high 
temperature to the ultracold regime. Molecules such as 
NH3, OH and NH have been cooled from room tempera- 
ture to the milliKelvin regime by a variety of methods in- 
cluding buffer-gas cooling [2(J, W$ an d Stark deceleration 
[22I [23. Directly cooled molecules have been success- 
fully trapped at temperatures around 10 mK, and there 
are a variety of proposals for ways to cool them further, 
including evaporative cooling, sympathetic cooling and 
cavity-assisted cooling 0, |25| . 

The possibility of controlling molecular interactions in 
the same way as atomic interactions is of great interest. 



Several groups have begun to explore the effects of ex- 
ternal fields on ultracold molecular collisions 26]. Volpi 
and Bohn [27| investigated collisions of 17 02 with He and 
found a very strong enhancement of spin-flipping cross 
sections even for weak magnetic fields. Krems et al. [2^] 
and Cybulski et al. |29| investigated spin-flipping colli- 
sions of NH with He and found a similar dependence on 
magnetic field. Krems and Dalgarno |30L |3l| elaborated 
the formal theory of scattering in a magnetic field for 
a variety of atom-molecule and molecule-molecule cases 
involving molecules in 2 S and 3 S states. Ticknor and 
Bohn 32] investigated OH+OH collisions and found that 
in this case magnetic fields could suppress inelastic col- 
lisions. Lara et al. [33l 13^ ] investigated the very com- 
plicated case of OH + Rb collisions, using basis sets de- 
signed to allow a magnetic field to be applied, though 
their initial calculations were for zero field. 

The effects of electric fields have also been investi- 
gated. Avdeenkov and Bohn [35l l3t| investigated OH- 
OH collisions in the presence of electric fields that caused 
alignment of the molecules. They identified novel field- 
linked states arising from long-range avoided crossings 
between effective potential curves in the presence of a 
field |M El S IS. Avdeenkov et al. gjj have also 
investigated the effects of very high electric fields on col- 
lisions of closed-shell molecules. Very recently, Tscher- 
bul and Krems ^l| have explored the effect of combined 
electric and magnetic fields on He-CaH collisions and ob- 
served significant suppression of spin-flipping transitions 
at high electric fields in the cold regime (~ 0.5 K). 

An important technique used to produce dimers in ul- 
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tracold atomic gases is magnetic tuning across Feshbach 
resonances 0, y| ■ A Feshbach resonance 0> OH occurs 
whenever a bound state associated with one potential 
curve lies above the threshold for another curve. The 
resonance thus corresponds to a level embedded in a con- 
tinuum, which is a quasibound state. In the atomic case, 
the thresholds that produce low-energy Feshbach reso- 
nances are associated with different hyperfine states of 
the interacting atoms. It is often possible to tune a reso- 
nance across threshold (from above or below) by applying 
a magnetic field. This produces an avoided crossing be- 
tween atomic and molecular states. If the magnetic field 
is tuned across the resonance slowly enough to follow the 
avoided crossing adiabatically, pairs of atoms can be con- 
verted into molecules or vice versa. 

Molecules have a much richer energy level structure 
than atoms, and there are many additional types of Fes- 
hbach resonance. In particular, rotational Feshbach res- 
onances can occur |44l 1 151 and have significant influences 
on ultracold molecular collisions 46] . Other small energy 
level splittings, such as spin-rotation and A-doubling, can 
also cause resonances in molecular scattering. It is of 
great interest to characterize such resonances and their 
field-dependence, both to understand their influence on 
collision cross sections and to prepare the ground for 
experiments that associate molecules by Feshbach res- 
onance tuning. 

The NH molecule is particularly topical in cold and ul- 
tracold molecule studies. It is a dipolar molecule with a 
3 £~ ground state, so it is both electrostatically and mag- 
netically trappable. It has been cooled by beam-loaded 
buffer-gas cooling |2l| and is a promising candidate for 
molecular beam deceleration and trapping |4jj. Krcms 
et al. |28| and Cybulski et al. j2(| have calculated poten- 
tial energy surfaces for He-NH and used them in scatter- 
ing calculations. Cybulski et al. also calculated zero- field 
bound states of the He-NH Van der Waals complex. An 
electronic excitation spectrum of the Van der Waals com- 
plex has been observed by Kerenskaya et al. • Soldan 
and Hutson 49] have investigated the interaction poten- 
tials for NH with Rb, and Dhont et al. 50] have devel- 
oped interaction potential energy surfaces for NH-NH in 
the singlet, triplet and quintet states. 

In the present article we describe the first calculations 
of the bound states of a Van der Waals complex in a mag- 
netic field. We show how such calculations can be used to 
locate zero-energy Feshbach resonances as a function of 
applied field. Elastic and inelastic cross sections can then 
be tuned by sweeping the field across the Feshbach res- 
onance. Such field sweeps could also be used to transfer 
unbound atom-molecule or molecule-molecule pairs into 
bound states of the corresponding complex. 

A remarkable conclusion of the present paper is that, in 
the presence of inelastic scattering, the poles in scattering 
lengths that characterize low-energy Feshbach resonances 
in atomic systems Q are dramatically suppressed and 
elastic and inelastic cross sections show relatively small 
peaks as resonances cross thresholds. The numerical re- 



sults obtained here allow us to test analytical formulae 
for this effect recently given by Hutson [5l| . 



II. METHODS FOR BOUND-STATE 
CALCULATIONS 

We consider the case of an NH molecule interacting 
with a He atom in the presence of a magnetic field. The 
Hamiltonian for this in Jacobi coordinates (R, 9) is 
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-H mon +H z +V(R,9), (1) 



where L 2 is the space-fixed operator for end-over-end ro- 
tation, -ffmon is the Hamiltonian for the NH monomer, Hz 
is the Zeeman interaction and V(R, 9) is the intermolecu- 
lar potential. For simplicity we consider the NH molecule 
to be a rigid rotor, but the generalization to include NH 
vibrations is straightforward. The NH monomer Hamil- 
tonian is therefore 

H mon = h- 2 b NU N 2 +H SN + H ss , (2) 

where &nh = 16.343 cm' 1 is the rotational constant of 
NH in its ground vibrational level [52| . 
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is the spin-rotation operator, and 
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is the spin-spin operator written in space-fixed coordi- 
nates [53. N and S are the operators for the rota- 
tional and spin angular momenta. The numerical val- 
ues for the spin-rotation and spin-spin constants are 
7 = -0.0055 cm" 1 and A ss = 0.920 cm" 1 

There are several basis sets that could be used to ex- 
pand the eigenfunctions of Eq. We consider two of 
them in the present work, which we refer to as the cou- 
pled and uncoupled basis sets. Both basis sets represent 
the end-over-end rotation with quantum numbers \LM£) , 
where L is a rotational quantum number and Ml is its 
projection onto the space-fixed Z axis. 

We use the convention that quantum numbers that de- 
scribe a monomer are represented with lower-case letters, 
and reserve capital letters to describe states of the com- 
plex as a whole. In the absence of a magnetic field (or 
a perturbing atom) , the rotational states of NH are ap- 
proximately described by quantum numbers n, s and j, 
where n represents the mechanical rotational of NH, s 
is the electron spin, and j is the vector sum of n and 
,s. In the coupled representation for the He-NH problem, 
we use basis functions \nsjmj)\LMi) that retain these 
monomer quantum numbers, with rrij the projection of j 
onto the space-fixed Z axis. In the uncoupled represen- 
tation, we use instead basis functions \nm n )\sm s )\LMi) , 
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where m n and m s are the projections of n and s individ- 
ually. 

In both basis sets we use, the matrix elements of L 2 
are diagonal in all quantum numbers and are simply 



h 2 L(L+l). The rotational part of the monomer Hamilto- 
nian is also diagonal, with matrix elements b^Yin(n + 1). 
The remaining matrix elements of the NH monomer 
Hamiltonian in the two basis sets are 



(nsjmjlHs^n'sj'm'j) = S nn ,6 jf 5 mjrn > 7 (-l) n+ ^+ s [n(n + l)(2n + l)s(s + l)(2s + 1)]* | * ™ { } ; (5) 
(nsjm^Hssln'sj'm'j) = 8 jf 6 mjm > ^A ss (-l) J +"' + " +s [(2n + l)(2n' + 1)]* ( J q o ){n s 2} ( 6 ) 

and 

(sm s I (nm„ | Hsn \ n'm' n ) \ sm' s ) 

+ (5 nn ,6 mnm , n±1 S msm , Tl ) 1 [n(n + 1) - <(r< ± 1)]* [s(s + 1) - m' s (m' s T 1)]* (7) 
(sm s \(nm n \H 3S \n'm' n )\sm' s ) = ^-X ss (-iy- m ^ \{2n + l)(2n' + 1)]* + l)(2a + 1)] f J § £ ) 

»{;i;}B-«-(- 



// 2 n' \ / s 2 s 
-m„ -g m' / I -m s g m' a 



(8) 



It may be noted that H mon is approximately diagonal in the coupled representation but not in the uncoupled repre- 
sentation. In the coupled representation, the only off-diagonal terms are matrix elements of Hgg that couple different 
rotational states with An = ±2. 

The Zeeman Hamiltonian for NH, neglecting rotational and anisotropic spin terms [54[ is 

Hz = . 9c /i B B • S, (9) 

where g c is the g-factor for the electron, /zb the Bohr magneton and B is the magnetic field vector. The matrix 
elements of this operator are 

(nsjmjlHzln'sj'm'j) = 6 nn ,S mjm ,g cf i B B(-l) n + s -^+ 1 [s(s + l)(2s + l)(2j + l)(2f + 1)]* 

3 1 3' \ h 3' n \ 



rrij ui j I 1 j s 1 
and 

(sm s \(nm n \Hz\n'rn' n )\sm' s ) = 5 nn tS mnm ' n 6 msm ' s g e fi B Bm a , (11) 

where the magnetic field direction has been chosen as the Z axis and B is the field strength. This is diagonal in the 
uncoupled representation but not in the coupled representation. 

The intermolecular potential is conveniently expanded in Legendre polynomials, 

V(R,9) = ^V A (i?)P A (cos0). (12) 

A 

The matrix elements of the Legendre polynomials in the coupled and uncoupled basis sets are 

{LM L \{nsj mj \Px{cos6)\n' sj'm'^L' M' L ) = [(2n+l)(2ra' + l)(2j + l)(2j' + 1)(2L+1)(2£' + 1)]* 

n X n' \ f L X L' \ v^/_ i^+i+i'+A+mA-Mi-mj 
^ l J 

L X L> \f j X f \ f j j> X 
—Ml —ni\ M' L J y —mj m\ m'j J \ n' n s 



(13) 
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This is off-diagonal in both representations. 

We solve the bound-state Hamiltonian by a coupled channel method (55|. Denoting a complete set of channel 
quantum numbers (n, s, j, TOj, L, M^) or (n, m n , s, m s , L, Ml) by i, we expand the total wavefunction 



(15) 



where (R,r) represents all coordinates except the intermolecular distance R and the channel functions <&i(R,f) are 
the corresponding basis functions of the coupled or uncoupled basis sets. Substituting this expansion into the total 
Schrodinger equation yields a set of coupled differential equations for the radial functions Xi {R) , 
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In the present work we solve the coupled equations to find 
bound states using the BOUND program |56||. which uses 
the algorithms described in ref . |5fJ For N channels there 
are N coupled equations. However, it is actually neces- 
sary to propagate a set of N linearly independent solu- 
tions, so x(R) is an NxN matrix. The log-derivative ma- 
trix Y — (dx/dR)x~ 1 is propagated outwards from i? m i n 
and inwards from i? max to a common matching point 
R m id in the classically allowed region of the potential us- 
ing Johnson's algorithm 57]. This is done for a series of 
trial energies. If the energy is an eigenvalue of the Hamil- 
tonian, the determinant of the log-derivative matching 
matrix Ymatch = Yin — Y out is zero. Bound states are lo- 
cated by searching for zeroes of eigenvalues of Ymatch as a 
function of energy. The log -derivative method provides a 
generalised node count |55l l58j which increases by one at 
each energy eigenvalue, and this allows us to use bisection 
to identify regions of energy that contain a bound state. 
The actual convergence on an energy eigenvalue uses the 
secant method, which gives quadratic convergence. 



Version 5 of the BOUND program [5(| contained an 
interface to allow new basis sets to be added, but could 
only handle basis sets in which the monomer Hamiltonian 
was diagonal. We have extended the program to remove 
this restriction and implemented the coupled and uncou- 
pled basis sets described above. We have also built in 
new loops over external fields to simplify calculations on 
Stark and Zeeman effects. 



III. RESULTS OF BOUND-STATE 
CALCULATIONS 

We have carried out bound-state calculations on 4 He- 
NH and 3 He-NH using the potential energy surface of 
Krems et al. (described in more detail as potential 2 
of Cybulski et al. [2{|). The basis set included all func- 
tions with n < 8 and L < 7. The coupled equations were 
propagated outwards from 1.8 A to 3.57 A and inwards 
from 16.0 A to 3.57 A using a log-derivative sector size 
of 0.025 A. WKB boundary conditions were applied in 
each channel at -R max to improve the convergence. 

In zero field, the levels of He-NH are characterized by 
the total angular momentum J7, which is the vector sum 
of j and L. The total parity is also conserved and is given 
by (— l)™ +i+1 . Since the lowest levels of NH have n = 0, 
j = s = 1 and J = L, L ± 1. The three different J levels 
corresponding to each value of L are very close together: 
the separation is only about 10~ 3 cm -1 for L = 1. The 
angular momentum coupling scheme corresponds to case 
(B) of Dubernet, Flower and Hutson |59|. 

The field-free energies of He-NH on this potential have 
been calculated previously by Cybulski et al. 29] . Our 
results agree with theirs to ±10~ 4 cm -1 for the levels 
with L — to 2, but are approximately 0.0014 cm -1 
lower for the L = 3 levels of 4 He-NH, which are bound 
by only 0.765 cm" 1 . We attribute the difference to lack 
of convergence of the radial basis set used in their calcu- 
lations. The results obtained from our program with the 
coupled and uncoupled basis sets are identical to ±10 -9 
cm -1 , which confirms the correctness of the code. 

In the presence of a field, J is very quickly destroyed. 
The only rigorously good quantum numbers in a mag- 
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FIG. 1: Bound-state energy levels for 3 He-NH for Mj = 
(upper panel) and Mj = — 1 (lower panel) as a function of 
magnetic field B. Levels of odd parity are shown in red and 
levels of even parity in green. Dissociation thresholds are 
shown as dotted blue lines and quasibound levels are shown 
as dashed lines. The arrows show the positions at which levels 
cross L = thresholds. 



netic field are parity and Mj = mj + Ml = fin + m s + 
Ml ■ The bound-state energies for levels correlating with 
n = for 3 He-NH with Mj = and —1 are shown in 
Fig. ^ Each level splits into components that can be la- 
belled with the approximate quantum numbers m n — 
and m s = 0, ±1. L remains an essentially good quantum 
number except in the vicinity of the avoided crossings. 
The major difference between M j values is that some 
m s levels are missing for Mj ^= 0, but there are also 
small shifts of all the energy levels. A similar plot for 
4 He-NH with Mj = is shown in Fig. H 

As expected, the Zeeman effect is generally quite linear 
for He-NH in the range of fields studied. This is because 
for n — the only off-diagonal terms are those of V (R, 9) 
and Hss that mix in excited n levels, and the spacing 
between the n = and n = 1 levels for NH is around 32.6 
cm -1 . Nevertheless, there are avoided crossings where 
levels of the same Mj and parity but different m s cross. 
Expanded views of the energy level diagrams in the region 



FIG. 2: Bound-state energy levels for 4 He-NH for Mj = 
as a function of magnetic field B. Levels of odd parity are 
shown in red and levels of even parity in green. Dissociation 
thresholds are shown as dotted blue lines and quasibound 
levels are shown as dashed lines. Avoided crossings are circled. 



of avoided crossings are shown for 4 He in Fig. [3] and for 
3 He in Fig. 0] ft may be seen that in this system the 
crossings are very tightly avoided, with spacings at the 
crossing points of AE = 5.2 x 1CU 4 cm -1 and AE = 
2.0 x 10" 3 cm" 1 for 4 He and AE = 1.5 x 1CT 3 cm" 1 for 
3 He. 

The reason that the crossings are so tightly avoided in 
this case is that there are no direct off-diagonal matrix 
elements between n = basis functions with different 
values of m s . The dominant coupling is a second-order 
one of the form 



(00lTO s |i? ss |2m rl lm' s )(2m n lm , s |T/ 2 |001m / s ; 



En 



E n 



(17) 



where basis functions are represented \nm n sm s ) and 
m n = m s — m' s . Such crossings will thus be more strongly 
avoided in systems with stronger anisotropy or smaller 
monomer rotational constant. 

If matrix elements of TJss off-diagonal in n are omitted, 
the avoided crossings still exist but are about a factor of 5 
tighter. Under these circumstances they are caused by a 
third-order mechanism in which potential couplings mix 
in basis functions with n > (but the same values of m s 
and m' s ) and these are connected by Hss- 

It should be noted that for levels with n > there are 
direct matrix elements of ifss and i?NS that connect lev- 
els with different m s but the same rrij . This will produce 
more strongly avoided crossings where the separation is 
simply proportional to Ass (or 7 for 2 S monomers). This 
will be particularly important for 16 02, which has an 
n = 1 ground state. 

The boundary conditions applied by the BOUND pro- 
gram are only correct for true bound states, below 
the lowest dissociation threshold. Above a threshold, 
the boundary conditions produce artificially quantized 
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FIG. 4: Avoided crossing between levels with m a = and 
for 3 He-NH as a function of magnetic field B. 



differences are that the energy E is above one or more 
thresholds Ei and that scattering boundary conditions 
are applied at long range. In the present work we solve 
the coupled equations for scattering using the MOLSCAT 
package [6(j. We have generalised MOLSCAT in the 
same way as BOUND to handle basis sets in which the 
monomer Hamiltonian is nondiagonal. 



A. Boundary conditions 



FIG. 3: Avoided crossing between levels with m a = +1 and 
m a = — 1 (upper panel) and between levels with m s = and 
m s = — I (lower panel) for 4 He-NH as a function of magnetic 
field B. 



states. These are of two types: states that are predomi- 
nantly in an open channel, which have no physical signifi- 
cance; and states that are predominantly in a closed chan- 
nel, which correspond closely to quasibound states of the 
real system but with only an approximate open-channel 
(dissociative) component. The two types are very easy 
to tell apart: the open-channel states are closely parallel 
to a lower threshold, and inspection of the wavefunction 
confirms their parentage. The open-channel states have 
been omitted from Figs. 1 and 2. The closed-channel 
states are shown as dashed lines above the lowest thresh- 
old. They allow us to estimate very precisely where a 
bound or quasibound state crosses a threshold, and thus 
where a zero-energy resonance is expected in the scatter- 
ing. This information will be used in section Ivl below. 



IV. METHODS FOR SCATTERING 
CALCULATIONS 

The coupled equations needed for scattering calcula- 
tions are identical to those for bound states. The only 



The usual procedure for obtaining the scattering ma- 
trix S from the log-derivative matrix Y has been de- 
scribed by Johnson |57| . Each channel is either open, 
E > Ei or closed, E < Ei. If matching takes place at a 
finite distance i? max where the wavefunction in the closed 
channels (E < Ei) has not decayed to zero, it is necessary 
to take account of the closed channels. The asymptotic 
form of the wavefunction is 

X(R) = J(R) + N(R)K, (18) 

where J(R) and N(R) are diagonal matrices made up 
Riccati-Bessel functions for open channels and modified 
spherical Bessel functions for closed channels. If there 
are N channels and N opcn open channels, the N x N 
log-derivative matrix is converted into an N x N real 
symmetric K matrix. The S matrix is then obtained from 
the open-open submatrix, K QO , using the relationship 

S = {I + iK 00 )- 1 {I-iK 00 ), (19) 

where I is an iVo pC n x N open unit matrix. However, the 
boundary conditions 1|18|) are appropriate only in a basis 
set in which both L 2 and the asymptotic Hamiltonian are 
diagonal. In our generalised version of MOLSCAT, the 
log-derivative matrix Y is propagated in the primitive 
basis set (which is nondiagonal) and then transformed at 
B max into a basis set that diagonalises H ulon + Hz ■ 
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For the simple case of He-NH, the eigenvalues of 
H m on + Hz are nondegenerate except at non-zero field. 
The transformation that diagonalises H mon + Hz is thus 
unique. However, in more complicated cases with two 
structured collision partners it will be necessary to trans- 
form to a basis set that diagonalises H m0 n,l + Hz,i and 
H m0 n,2 + Hz.2 separately. Additional degeneracies arise 
at zero field, and resolving them may require diagonali- 
sation of another operator such as jz ■ 

B. Cross sections 



However, this relies on a Taylor series expansion of 
exp(2i(5) that is valid only when 5 <C 1. Across an elastic 
scattering resonance, 6 changes by ir even at limitingly 
low energy, so Eq. 1261 is inappropriate. In the present 
work we obtain scattering lengths numerically either by 
converting the low-energy S-matrix elements to complex 
phases using Eq. |^ and then the definition l|22|) or by 
using the equivalent identity 



The cross section for a transition i — ► / from initial 
state i to final state / is obtained from the square of the 
corresponding T matrix element, 



(20) 



where k is the incoming wave vector, k 2 — 2fx{E— Ei)/Ti 2 , 
and Tif = Sif — Sif. In general it is necessary to sum 
over all channels corresponding to the monomer levels of 
interest and over S matrices obtained for different values 
of Mj and parity. 



C. Scattering lengths 

In the ultracold regime, scattering properties are often 
described in terms of a complex scattering length a = 
a — \(3 0, ■ The diagonal S-matrix element in the 
incoming channel may be written in terms of a complex 
phase shift 5 |63| . 



S o = exp(2i(5) 
and the complex scattering length is defined by 

— tan (5 



Equivalently, 



Sn 



ike 



1 + ika 



(21) 



(22) 



(23) 



The scattering length becomes constant at limitingly low 
energy. The elastic and total inelastic cross sections are 
M 



47r|a| : 



Cel = 



and 



Cinel 



1 + k 2 \a\ 2 + 2k(3 

47T/3 



k(l + k 2 \a\ 2 + 2k(3)' 

The scattering length is often given as 

,. 1 T 0Q 
a = urn — — — . 

fe^o 2i k 



(24) 



(25) 



(26) 



D. Resonant behavior 

If there is only one open channel, then the phase shift 
5 is real and its behavior is sufficient to characterize a 
resonance. It follows a Breit-Wigner form as a function 
of energy 



8(E) 



:>b g 



tan 



r, 



2[E Ies — E) 



(28) 



where <5b g is a slowly varying background term, E les is the 
resonance position and Te is its width (in energy space) . 
This corresponds to the S matrix element describing a 
circle of radius 1 in the complex plane. In general the 
parameters 5bg, -E*es and Te are slow functions of en- 
ergy, but this is neglected in the present work apart from 
threshold behaviour. 

The resonance position and the threshold energy 
-E"thresh are both functions of magnetic field, 



dE r . 



dB 



Mr 



and 



dE, 



thresh 



dB 



Mthrcsh ■ 



(29) 



We define B res (E) as the field at which E res — E. Fol- 
low kinetic energies, the width Te also depends on the 
field through its threshold dependence [65|. 



F E (E kin ) = 2kj E , 



(30) 



with constant reduced width *y E . As a function of mag- 
netic field at constant kinetic energy, the phase shift thus 
follows a form similar to Eq. 



6(B) = dbg + tan~ 



Fs^kin) 



2{B rcs {E) - B) 



(31) 



The width TsiE^) is a signed quantity that is negative 
if the bound state tunes upwards through the energy of 
interest and positive if it tunes downwards, 



Ts(i?kin) 



kin) 



Mthrcsh Mrcs 



(32) 



The background phase shift (5b g goes to zero as k — > 
according to Eq. 1221 (with constant finite a^g), but the 
resonant term still exists. The scattering length passes 
through a pole when <5 = (n + h) ft- The scattering 
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length follows the formula commonly used in atomic scat- 
tering 66], 



a(B) = a hg 



1 - 



B - B res (E) 



(33) 



where Ab = — rs/2fcab g = — Is/o-bg and is independent 
of k near threshold. 

When there are several open channels, 5 is in gen- 
eral complex. The quantity that then follows the Brcit- 
Wigner form (|28|l or (|31|1 is the S-matrix eigenphase sum 
|44j , which is the sum of phases of the eigenvalues of the 
S matrix. The eigenphase sum is real, because the S ma- 
trix is unitary, so that all its eigenvalues have modulus 1. 
In practice the eigenphase sum is most conveniently cal- 
culated by diagonalising the real symmetric matrix K OQ 
and summing the inverse tangents of its eigenvalues. 

If there is more than one open channel, the individual 
S matrix elements still describe circles in the complex 
plane, 



bg.ii' 



^9Ei9Ei 



E-E rcs +iT E /2' 



(34) 



where g E i is complex. The partial width for channel i 
is T E i — \gEi\ 2 and the radius of the circle in Su> is 
\9EigEi' |/r_E- The analogous expression as a function of 
magnetic field at constant kinetic energy is 



Sut — S 



bg,ii' 



igBigBi 



B-B Tes +iT B /2> 



(35) 



where the energy-dependence of B les and T E has been 
omitted to simplify notation. If Teo < ^e (or equiv- 
alently \Tbo\ < \^b\), the scattering length does not 
pass through a pole. In the low-energy threshold regime, 
Tbo is proportional to k and we may define an energy- 
independent reduced width jbo, 



Tbo = ^iBok. 



(36) 



However, Tb has inelastic contributions T^ el that are 
essentially energy- independent. Hutson has defined 
a resonant scattering length, 



27. 



BO 



1 B 



(37) 



that characterizes the strength of the resonant contribu- 
tion to the scattering at low energy. If ka rcs <C 1, Sqq 
describes a circle of radius 2fca ros in the complex plane 
as a resonance is tuned through threshold and the real 
part of the scattering length oscillates by ±a res /2. 

As will be seen below, for He + NH the inelastic 
scattering strongly suppresses the pole in the scatter- 
ing length and the resonant oscillation in the scattering 
length is of quite small amplitude. The corresponding 
oscillations in the elastic and inelastic cross sections are 
also relatively weak. 



V. RESULTS OF SCATTERING 
CALCULATIONS 

We have carried out scattering calculations on 4 He-NH 
and 3 He-NH using the potential energy surface of Krems 
et al. 28}. These calculations used a reduced basis set 
of functions with n < 4 and L < 5 to allow comparison 
with previous studies. The coupled equations were prop- 
agated outwards from 1.7 A to 120.0 A using Johnson's 
log-derivative algorithm [57j with a sector size of 0.025 

A. 

We have verified that the new code gives identical scat- 
tering results for He-NH with the coupled and uncoupled 
basis sets. This confirms the correctness of the coding 
for the basis sets and the extraction of S matrices. Our 
program also gives identical results to refs.l28l andl29l for 
scattering in a magnetic field. 

It is of great interest to investigate the effect of zero- 
energy Feshbach resonances on low-energy molecular 
scattering as a function of magnetic field. However, a 
major problem is that the resonances can be very narrow, 
and locating the fields at which they occur is difficult and 
time-consuming. For example, in the He-NH problem 
we expect the coupling between bound and continuum 
states with different mj to be comparable to or smaller 
than that between bound states of different mj. To a 
first approximation, the latter is the energy separation of 
the bound-state avoided crossings, which is around 10 -3 
cm . Since the energy of a state with m s — ±1 tunes by 
about 10~ 4 cm _1 /G, we expect the Feshbach resonances 
to be less (perhaps much less) than 10 G wide. 

Fortunately, the bound-state capability described 
above provides a solution to this problem: we can ex- 
trapolate the bound-state energies to calculate the field 
at which they cross a threshold and then scan across a 
small range of fields in the vicinity. As an example, we 
searched for Feshbach resonances in 3 He-NH. It may be 
seen in Fig. 2] that there is a 3 He-NH bound state with 
m s = +1 that crosses the m s = — 1 threshold at about 
7200 G and the m s = threshold at about 14300 G. 
Careful extrapolation of the bound-state energies with 
the basis set used for the scattering calculations gives 
more precise field estimates of 7168.750 G (for Mj = — 1) 
and 14340.36 G (for Mj = 0) respectively. 

Care must be taken to use bound-state calculations 
for the correct values of Mj and parity. In the present 
case, we want s-wave resonances, with L = in the in- 
coming channel. This requires Mj = rrij and even par- 
ity. Bound-state calculations with different values of Mj 
produce energies that cross threshold at different values 
of the field. In He-NH, which is a very weakly coupled 
system with almost linear Zeeman effects, the fields are 
only slightly different; for example, the crossing with the 
m s = threshold occurs at 14345.40 G for Mj = -1. 
Nevertheless, the difference can easily be enough to miss 
the resonance, and in more strongly coupled systems will 
be crucial. 

For 3 He colliding with NH (m s = -1) near 7169 G, 
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7168.750 7168.755 7168.760 

Magnetic Field (G) 



7168.765 



3e+04 



2e+04 

■< 



a 1e+04 - 



0e+00 



-1e+04 - 



-2e+04 - 



-3e+04 



. . . I 

E kin = 1.0 

J 




1 


- ■ ' ' ' 



7168.7550 



7168.7555 
Magnetic Field (G) 



7168.7560 



FIG. 5: Elastic (m, = -1 — > -1) cross section for He-NH 
collisions in the vicinity of an elastic Feshbach resonance at a 
kinetic energy of 10 -6 K. 



FIG. 7: Scattering length for elastic He-NH collisions in the 
vicinity of an elastic Feshbach resonance, from calculations at 
a kinetic energy of 10 ~ 6 K. 



9.0 



1 ■ 1 ■ 1 















I.I.I 



7168.752 7168.754 7168.756 7168.758 
Magnetic Field (G) 

FIG. 6: Phase shift for elastic 3 He-NH collisions in the vicinity 
of an elastic Feshbach resonance at a kinetic energy of 10 -6 K. 



only elastic scattering can occur. For Mj = —1 and 
even parity, the present basis set gives 3 open channels 
with L = 0, 2 and 4. Scattering into the L > chan- 
nels is strongly suppressed by the centrifugal barriers. 
The elastic cross section is shown as a function of field in 
Fig. [S] and the corresponding eigenphase sum and scat- 
tering length are shown in Figs. and The peak in 
the cross section for kinetic energy E^ n = 1CP 6 K is 
close to the value of 4ir/k 2 = 1.2 x 10 8 A 2 characteris- 
tic of a pole in the scattering length. The peak shows 
an asymmetric Fano lineshape , with interference be- 
tween a background term and a resonant term that inter- 
fere constructively on the low-field side of the resonance 
and destructively on the high-field side. The diagonal S 
matrix element for L = describes a circle of radius 1 
in the complex plane as the field is ramped across the 
resonance, as shown in Fig. [S] Fitting the eigenphase 
sum to Eq. E2 gives 5 bg < 0.001, B Tes = 7168.7555 G 
and Tb = —1.65 x 10~ 5 G, while fitting the scattering 
length to Eq.Elgives a bg = 3.19 A, B rcs = 7168.7555 G 
and A B = 8.04 x 1CT 3 G. The widths r B and A B are 




-1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 



Re [S 00 ] 

FIG. 8: The circle of radius 1 in the complex plane described 
by the L — diagonal S matrix element for elastic scattering 
in the vicinity of an elastic Feshbach resonance at a kinetic 
energy of 10 -6 K. The cross shows the value far from reso- 
nance. 

related by Tb/2 = —ka^ g As- It may be seen that the 
bound-state calculation does indeed give a very precise 
estimate of the position of the zero-energy resonance. 

For 3 He colliding with NH (m s = 0) near 14300 G, 
both elastic and inelastic scattering are possible. For 
Mj = and even parity, our basis set gives 3 elastic 
channels (m s = with L = 0, 2 and 4) and 2 inelastic 
channels (m s = — 1 with L = 2 and 4). However, the 
elastic channels with L > make no significant contri- 
butions at ultralow energies. Fig. [5] shows scans of the 
elastic (m s =0^0) and total inelastic (m s = — ► — 1) 
cross sections for kinetic energies of 10 ~ 6 K and 10~ 3 K. 
Once again the bound-state calculation gives a very pre- 
cise estimate of the position of the zero-energy resonance. 
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E kin = 1.0 nK 
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m„ = 0, m ' = -1 



14340.2 14340.3 14340.4 

Magnetic Field (G) 



1e+00 
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14347.7 14347.8 14347.9 

Magnetic Field (G) 



14340.5 




14348.0 



FIG. 9: Elastic (m s = — » 0) and inelastic (m s = — > 
— 1) cross sections for 3 He-NH collisions in the vicinity of an 
inelastic Feshbach resonance at kinetic energies of 10 -6 K 
(upper panel) and 10 -3 K (lower panel). The lines show the 
results of Eqs. 1591 and 1301 




1.30 14340.34 14340.38 

Magnetic Field (G) 



14340.34 14340.38 
Magnetic Field (G) 



14340.42 




14340.42 



FIG. 10: Real and imaginary parts of the scattering length 
for 3 He-NH collisions in the vicinity of an inelastic Feshbach 
resonance, from calculations at a kinetic energy of 10~ 6 K. 
Upper panel: real part; lower panel: imaginary part. The 
lines show the results of Eq. 1381 



At 10 -3 K the resonance is shifted slightly because a dif- 
ferent field is needed to bring the bound state into reso- 
nance with the larger total energy. Apart from the shift, 
however, the cross sections behave as expected from the 
Wigner threshold laws: the elastic cross section is almost 
unchanged and the inelastic cross section scales with k^ 1 . 

A notable feature of Fig. [5] is that the elastic cross 
section does not peak at the very high value of An/k 2 
characteristic of a pole in the scattering length. The real 
and imaginary parts of the scattering length are shown 
in Fig.^3 Instead of rising to oo, the scattering length 
oscillates with a peak at less than +8 A. The eigenphase 
sum shows a sharp drop through ir as shown in Fig. 1111 
and fitting to Eg. I3T1 rives B les = 14340.371 G and T B = 
— 5.72x 10~ 3 G. However, the phase change is distributed 
between several diagonal matrix elements. 

The S matrix elements are shown explicitly for two 
different energies, Skin = 10~ 6 K and 4 x 10~ 6 K, in 
Figs. 1121 and 1131 The corresponding values of k differ by 
a factor of 2. It may be seen in Fig. ^| that the elastic S 
matrix element for the low-energy incoming channel de- 
scribes a small circle with a radius that depends linearly 
on k. For £"kin = 10~ 6 K, the radius is only about 0.003. 



The other diagonal elements describe much larger circles 
that are almost independent of k, with background (non- 
resonant) S matrix elements (shown by crosses) that are 
far from 1: these channels have substantial kinetic energy 
so are not governed by the Wigner threshold laws. The 
inelastic S matrix elements are shown in Fig. 1131 and it 
may be seen that those for incoming channel describe 
circles with radius proportional to k 1 ^ 2 . 

Since Sqo describes a small circle in the complex plane, 
the corresponding complex phase S given by Ea. l21l shows 
a relatively small oscillation that does not pass through 
±7t/2. There is thus an oscillation but no pole in the scat- 
tering length a (and no discontinuity in its sign). This 
explains why the peak in the elastic scattering cross sec- 
tion is much lower than the value of 4ir/k 2 expected when 
there is a pole in a. 

The behavior observed here may be quantified in terms 
of the theory of Hutson [5l|. Fitting the individual S 
matrix elements at i^in = 10~ 6 K to Eq. [35] gives a 
partial width Tbo = —1-646 x 10~ 5 G for the incoming 
channel and r fll = -5.693 x 10" 3 G and T B2 = -1.308 x 
10~ 5 G for the two inelastic channels. We have verified 
numerically that T^o is proportional to k and Tbi and 
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TABLE I: Positions (A) and heights (K) of centrifugal barri 



9.0 ■ 



14340.30 14340.34 14340.38 
Magnetic Field (G) 



14340.42 



FIG. 11: Eigenphase sum for inelastic He-NH collisions in 
the vicinity of an inelastic Feshbach resonance at a kinetic 
energy of 10 -6 K. 



T B 2 are independent of it. With k = 3.219 x 1(T 4 A" 1 , 
this gives 7bo = 0.02557 G A and a rcs = 8.96 A. The 
peaks in the scattering lengths in Fig. 1 101 are at a ma x = 
7.67 A, a min = -1.28 A, and /3 max = 8.96 A, which 
correspond to the theoretical predictions jHJ], a™f* = 
ct bg ± a ros /2 and /3 max = a rcs with a bg = 3.19 A. 

The shapes of the peaks in scattering lengths and cross 
sections also correspond very accurately to the theoreti- 
cal profiles [5l"|. 





3 


He-NH 


4 


He-NH 


L 




V(-Rmax) 


-Rmax 


V(R W ) 


1 


9.6 


0.14 


10.2 


0.10 


2 


7.3 


0.73 


7.7 


0.52 


3 


6.1 


2.06 


6.5 


1.46 


4 


5.4 


4.43 


5.7 


3.14 



Table [fl Because of this, no zero-energy resonance was 
observed for Mj = —1 near 14345 G, even though there 
is a state that crosses the threshold there. 

We should emphasize the advantage of using coupled- 
channel methods rather than basis-set methods for bound 
states when attempting to locate zero-energy Feshbach 
resonances. Basis-set methods lose accuracy close to dis- 
sociation because of the difficulty of representing near- 
dissociation and continuum functions with basis sets. In 
coupled-channel methods, by contrast, the behavior of 
the wavefunction at long range can be built in by apply- 
ing WKB boundary conditions at i? max - Because of this, 
the approximations in the bound-state calculations are 
very similar to those in scattering calculations with the 
same basis set, and as we have shown this allows very 
precise estimates of resonance positions. 



a(B) 

<Jel{B) 



Obg 



2(B-5 rcs )/r£ cl + i' 



4n\a\ = 47r 



Che 



(38) 

2 



2{B - B ICS )/r^ ei + i 



47T/3 _ 47T Ores 

k ' k [2(B-B ICS )/r^] 2 + l 



(39) 
(40) 



The resulting profiles are shown in Figs.l9land ll(JI using 
the parameters determined from the S matrix elements 
without refitting to the scattering lengths and cross sec- 
tions. 

It is remarkable that, even in a system as weakly 
coupled as He-NH, the behavior of the cross sections 
and scattering lengths is so different from the elastic 
case. It has been common in theoretical studies of atom- 
atom scattering to model Feshbach resonances using a 2- 
channel treatment with one closed and one open channel. 
The present results make it clear that such an approxi- 
mation may miss essential features of the physics, and in 
particular may predict unphysical poles in the scattering 
length. 

Feshbach resonances with L > in the incoming chan- 
nel can also occur, but are suppressed at ultralow energies 
because of the centrifugal barrier. 3 He-NH has a very low 
reduced mass, and with a long-range potential given by 
—CqR~ 6 with Cq = 7.84_E/j <2q the heights and positions 
of the centrifugal barriers for L = 1 to 4 are given in 



VI. RELATIVE MERITS OF DIFFERENT BASIS 

SETS 



The coupled and uncoupled basis sets give identical 
results for bound-state energies and for collision proper- 
ties, and in this sense they are equivalent. However, the 
uncoupled basis set is a little easier to program and its 
matrix elements are easier to generalize to more compli- 
cated cases (such as those involving nuclear spin or two 
structured monomers). In addition, the uncoupled basis 
set gives a much simpler representation of the wavefunc- 
tions for bound states and resonances in any significant 
magnetic field. We therefore intend to use uncoupled 
basis sets in future work. 



VII. CONCLUSIONS 

We have modified the BOUND and MOLSCAT pack- 
ages to allow the use of basis sets in which the asymptotic 
Hamiltonian is non-diagonal and used the new capabil- 
ity to perform bound-state and scattering calculations on 
He-NH in the presence of a magnetic field. The bound- 
state capability makes it possible to locate zero-energy 
Feshbach resonances as a function of magnetic field even 
when they are very narrow. The new capability pro- 
vides a very straightforward way to program new cou- 
pling cases and collision types, and in future work we 



12 



will use it to investigate low-energy scattering of systems 
containing two structured monomers in electric and mag- 
netic fields. 

For He-NH, we have located two zero-energy Feshbach 
resonances involving an m s = +1 level of the He-NH 
complex tuned through the m s = — 1 and m s — thresh- 
olds. The two resonances show very different behavior 
as a function of magnetic field. For the resonance at 
the m s = — 1 threshold, only elastic scattering is possi- 
ble and the resonance shows the classic behavior famil- 
iar from ultracold atom-atom scattering. The scattering 
length passes through a pole at resonance and the elastic 
cross section shows a very large peak. For the resonance 
at the m s — threshold, however, inelastic scattering is 
also possible. The pole in the scattering length is dra- 
matically suppressed, and the peak in the elastic cross 
section is far smaller than would be expected if a pole 
were present. 

Our results provide a numerical demonstration of the 
effects recently predicted by Hutson jHJ, who param- 
eterized the strength of the resonant contribution with 
a resonant scattering length a rcs . For the inelastic reso- 
nance that we have characterized in He-NH, the resonant 
scattering length is only 8.96 A. 

The suppression of resonant peaks in cross sections by 



inelastic scattering will be a very general effect. Most 
molecular systems have stronger inelasticity than He-NH, 
and will have even smaller resonant scattering lengths. 
In such cases the peaks will be even more strongly sup- 
pressed than here. The effect explains previously puz- 
zling results in Na + Na2 H^] and Li + Li2 |0 scattering, 
where it was found that cross sections for ultracold colli- 
sions of vibrationally excited molecules were only weakly 
dependent on potential parameters and showed no sharp 
peaks when Feshbach resonances were tuned through zero 
energy. The lack of such peaks is now seen to be due to 
suppression by inelastic processes. This suggests that 
calculations on ultracold collisions of molecules in the 
presence of inelasticity may be much less sensitive to po- 
tential details than was previously expected. 
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FIG. 12: The circles described by diagonal S matrix ele- 
ments in the presence of inelastic scattering at Ey^n = 1CP 6 K 
(green) and 4 x 10~ 6 K (red). Upper panel: incoming wave 
(channel 0: m a = 0, L = 0); center panel: (channel 1: 
m s — —1, L = 2); lower panel: (channel 2: m s — —1, L = 4). 
The crosses show values far from resonance. 
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FIG. 13: The circles described by the off-diagonal S matrix 
elements for 3-channel scattering at -Ekm = 1CP 6 K (green) 
and 4 x 10~ 6 K (red). Channels are labeled as in Fig.lDl The 
crosses show values far from resonance. 



